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The  simulation  results  of  a  one-dimensional  (ID)  direct  methanol  fuel  cell  (DMFC)  model  are  compared 
with  the  current  density  and  methanol-crossover  data  that  are  experimentally  measured  under  several 
different  cell  designs  and  operating  conditions.  No  fitting  parameters  are  employed  for  the  compari¬ 
son  and  model  input  parameters  obtained  from  the  literature  are  consistently  used  for  all  the  cases  of 
comparison.  The  numerical  predictions  agree  well  with  the  experimental  data  and  the  ID  DMFC  model 
successfully  captures  key  experimental  trends  that  are  observed  in  the  cell  current  density  and  methanol- 
crossover  data.  This  clearly  illustrates  that  the  present  DMFC  model  can  be  applicable  for  optimizing  DMFC 
component  designs  and  operating  conditions.  In  addition,  the  model  simulations  further  indicate  that  the 
reduction  of  the  methanol  concentration  in  the  anode  catalyst  layer  is  critical  to  simultaneously  suppress 
both  the  electro-osmotic  drag  (EOD)  and  the  diffusion  aspects  of  methanol  crossover. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Within  the  last  decade,  portable  electronic  devices  such  as 
mobile  phones,  portable  cameras,  and  notebook  computers  have 
become  increasingly  miniaturized  and  multi-functionalized,  which 
generates  an  enormous  demand  for  superior  capacity  recharge¬ 
able  batteries  to  power  these  devices.  As  a  result,  the  use  of  these 
mobile  electronic  devices  has  been  limited  by  battery  performance 
that  suffers  from  short  operation  times  and  long  recharge  times. 
In  this  context,  direct  methanol  fuel  cells  (DMFCs)  using  polymer 
electrolyte  membranes  have  received  considerable  attention  for 
replacing  batteries  for  portable  electronics  devices  due  to  their 
potentially  higher  energy  specific  and  nearly  zero  recharge  time. 

Although  substantial  efforts  have  been  made  to  advance  DMFC 
technologies,  there  still  exist  several  issues  that  currently  prevent 
the  appearance  of  DMFCs  in  the  market.  These  are  related  to  the 
low  activity  of  electrocatalysts  for  the  methanol  oxidation  reac¬ 
tion  (MOR)  at  the  anode,  substantial  methanol  crossover  from  the 
anode  to  the  cathode,  and  water  imbalance  between  the  anode  and 
the  cathode  (the  requirement  for  anode  water  and  the  flooding  of 
cathode  water).  Furthermore,  to  exceed  the  energy  specific  of  Li-ion 
batteries,  it  is  necessary  to  feed  highly  concentrated  methanol  fuel 
directly  into  a  DMFC  stack;  this  worsens  the  methanol  crossover 
problem  and  renders  water  management  in  DMFCs  even  more  dif¬ 
ficult  [1]. 
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An  implication  of  the  above  observations  is  that  DFMC  com¬ 
ponent  design,  in  terms  of  physical  dimensions,  materials,  etc., 
is  extremely  important  to  suppress  methanol  crossover  from  the 
anode  to  the  cathode.  Therefore,  many  researchers  have  attempted 
to  measure  experimentally  the  degree  of  methanol  crossover  under 
various  DMFC  conditions  and  its  effects  on  DMFC  performance. 
Heinzel  and  Barragan  [2]  comprehensively  reviewed  the  main  fea¬ 
tures  of  methanol  transfer  through  Nation®  and  other  types  of 
membranes  in  DMFCs.  They  summarized  the  influence  of  various 
DMFC  designs  and  operating  factors  such  as  methanol  concen¬ 
tration,  pressure,  temperature,  membrane  thickness  and  catalyst 
morphology  on  methanol  crossover  and  the  resulting  cell  perfor¬ 
mance.  Ren  et  al.  [3  ]  studied  the  effects  of  methanol  crossover  under 
different  membrane  thicknesses  and  temperatures.  They  mea¬ 
sured  methanol  permeation  rates  across  Nafion®  112,115  and  117 
membranes  and  found  that  the  methanol  flux  rate  increases  with 
temperature.  Hikita  et  al.  [4]  estimated  the  amount  of  methanol 
crossover  by  measuring  the  species  composition  in  the  cathode 
exhaust  gas  under  different  methanol  feed  concentrations  (3,  6 
and  9  vol.  %)  and  membrane  thicknesses  (Nafion®  112,  115  and 
117).  They  found  that  the  use  of  a  thicker  membrane  successfully 
lowers  the  methanol  crossover  flux  but  it  also  lowers  the  cell  perfor¬ 
mance  due  to  the  larger  ohmic  resistance.  Ravikumar  and  Shukla  [  5  ] 
demonstrated  that  the  cathode  polarization  loss  increases  consid¬ 
erably  with  cell  temperature  and  methanol  feed  concentration  due 
to  the  enhanced  methanol  crossover.  Wang  et  al.  [6]  experimentally 
estimated  the  methanol  crossover  rates  and  found  that  methanol 
crossover  results  in  a  mixed  potential  and  poisoning  of  Pt  catalysts 
that  significantly  increases  the  polarization  loss  on  the  cathode. 
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Nomenclature 

c  molar  concentration,  mol  m-3 

D  diffusivity,  m2  s-1 

F  Faraday’s  constant,  96487  C  mol-1 

/  current  density,  mAcm-1 

i reference  exchange  current  density,  mA  cm-1 
I<  absolute  permeability  of  porous  media,  m2 

I<mem  membrane  hydraulic  permeability,  m2 
I<reac  reaction  rate,  s-1 

M  molar  mass,  kg  mol-1 

N  molar  flux,  mol  m-2  s-1 

nd  electro-osmotic  drag  coefficient 

P  pressure,  Pa 

R  universal  gas  constant,  8.314  J  mol-1  K-1 

r  pore  radii,  m 

Rcont  contact  resistance,  m£2  cm-1 

s  liquid  water  saturation 

T  operating  temperature,  K 

U°  standard  equilibrium  potential,  V 

Vceu  cell  voltage,  V 

Shortcuts 

aCL  anode  catalyst  layer 

BL  backing  layer 

cCL  cathode  catalyst  layer 

CCM  catalyst-coated  membrane 

MEA  membrane  electrode  assembly 

MEM  membrane 

Superscripts 

cons  consumption 

diff  diffusion 

eff  effective 

EOD  electro-osmotic  drag 

8  gas 

l  liquid 

prod  production 

Subscripts 
a  anode 

BC  boundary  condition 

c  cathode 

CO2  carbon  dioxide 

in  inlet 

k  species  phase 

mem  membrane 

M  methanol 

02  oxygen 

ref  reference 

t  total 

W  water 

xover  crossover 

Greek  symbols 

y  reaction  symmetry 

8  thickness,  m 

s  porosity 

r\  efficiency 

6  contact  angle,0 

k  proton  conductivity 

A  water  content 

ji  viscosity,  kg  cm-1  s-1 

£  stoichiometric  ratio 


p  density,  kg  m-3 

a  surface  tension,  N  m-1 

r  tortuosity 


As  mentioned  above,  it  has  been  experimentally  proven  that 
methanol  crossover  is  a  key  issue  in  DMFCs  and  that  there  are 
a  number  of  design  and  operating  parameters  that  influence  the 
degree  of  methanol  crossover  and  the  resulting  cell  performance. 
That  indicates  that  a  DMFC  demands  a  high  degree  of  optimization. 
For  the  successful  design  and  operation  of  a  DMFC  system,  a  better 
understanding  of  methanol-crossover  phenomena  is  essential  and 
can  most  likely  be  achieved  by  theoretical  DMFC  modeling,  as  well 
as  experimental  investigations. 

A  number  of  DMFC  models  have  been  developed  for  the  purpose 
of  obtaining  a  basic  understanding  of  key  physical  phenomena  in 
DMFCs  and  optimizing  cell  designs  and  operating  conditions.  These 
DMFC  models  range  from  simple  one-dimensional  mathematical 
models  [7-12]  to  complex  multi-dimensional,  CFD  models  [13-18]. 
Despite  the  fact  that  rapid  advances  in  DMFC  modeling  in  the  last 
decade  have  yielded  substantial  predictive  capabilities,  efforts  for 
model  validation  against  experimental  data  have  been  relatively 
limited  to  quite  a  narrow  range  of  DMFC  operating  conditions. 

Scott  et  al.  [7]  developed  a  single-phase,  semi-empirical  ID 
DMFC  model  and  investigated  the  impact  of  methanol  crossover 
on  cell  performance.  In  their  study,  the  results  of  the  simula¬ 
tion  model  were  compared  with  experimental  data  that  were 
measured  under  the  operating  conditions  of  a  methanol  feed  con¬ 
centration  of  1  M  at  the  anode  and  three  values  of  the  oxygen 
pressure  on  the  cathode  (1 , 2  and  3  bars).  They  reported  good  agree¬ 
ment  between  the  model  calculations  and  experimental  data  only 
over  a  range  of  moderate  current  densities  (350-400 mAcm-2); 
considerable  deviations  occur  at  low  and  high  current  densi¬ 
ties.  Murgia  et  al.  [8]  developed  a  one-dimensional,  two-phase 
model  for  a  liquid-feed  DMFC  and  intensively  studied  methanol- 
crossover  phenomena  with  different  methanol  feed  concentrations 
and  membranes.  They  successfully  validated  their  model  against 
the  experimental  data  measured  by  Hikita  et  al.  [4]  under  various 
methanol  feed  concentrations  (0.5,  1.0,  1.5  M).  In  the  model  vali¬ 
dation,  however,  several  key  parameters  related  to  water  flooding 
and  condensation/evaporation  were  considered  to  be  adjustable 
parameters  for  a  better  fit  with  the  experimental  data.  Chen  and 
Zhao  [9]  developed  a  ID,  non-isothermal  DMFC  model  and  investi¬ 
gated  thermal  effects  on  methanol  crossover  and  cell  performance. 
They  did  not,  however,  attempt  to  experimentally  validate  their 
model.  Recently,  Kareemulla  and  Jayanti  [10]  introduced  a  ID, 
single-phase  DMFC  model  and  studied  the  limiting  current  density 
and  methanol-crossover  phenomena  in  DMFCs.  They  decomposed 
the  methanol  crossover  flux  into  electro-osmotic  and  diffusion- 
driven  methanol  fluxes  and  intensively  analyzed  the  behaviour 
of  methanol  crossover  over  a  wide  range  of  methanol  feed  con¬ 
centrations  (0.125-5.0  M).  Their  model  validation,  however,  was 
conducted  under  a  narrow  range  of  methanol  feed  concentra¬ 
tions,  i.e.,  0.125,  0.25  and  0.5  M.  Oliveira  et  al.  [11]  presented  a 
one-dimensional,  single-phase  DMFC  model  and  precisely  stud¬ 
ied  methanol  crossover  and  water  transport  across  the  membrane. 
They  compared  simulation  results  with  experimental  data  that 
were  measured  at  low  methanol  feed  concentrations  (0.1,  0.2,  and 
0.5  M)  and  a  fixed  temperature,  viz.,  70  °C.  Eccarius  et  al.  [12]  for¬ 
mulated  a  one-dimensional  methanol  crossover  model  to  describe 
methanol  crossover  from  the  anode  to  the  cathode.  The  model 
calculation  results  were  compared  with  methanol  crossover  cur¬ 
rent  density  data  that  were  measured  under  various  methanol 
feed  concentrations,  temperatures,  and  membrane  thicknesses. 
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In  their  work,  however,  the  electro-osmotic  drag  (EOD)  coeffi¬ 
cient  and  the  methanol  concentration  in  the  cathode  catalyst 
layer  were  adjustable  parameters  for  fitting  with  the  experimental 
data. 

Multi-dimensional  DMFC  models  have  been  presented  by  many 
DMFC  modeling  groups.  Wang  and  Wang  [13]  presented  a  2D, 
two-phase  mixture  model  for  a  liquid-feed  DMFC.  The  model  pre¬ 
dictions  were  compared  with  experimental  data  that  were  obtained 
from  a  5  cm2  cell  under  various  operating  conditions  (1  and  2M 
methanol  concentrations;  50  and  80  °C  operating  temperatures; 
different  anode  and  cathode  flow  rates).  Divisek  et  al.  [14]  devel¬ 
oped  a  2D,  two-phase  DMFC  model  in  which  the  effects  of  both 
oxygen  and  methanol  electrochemical  kinetics  and  methanol  per¬ 
meation  across  the  membrane  were  considered.  Their  model  was 
experimentally  validated  under  various  methanol  feed  concentra¬ 
tions  (0.5,  1  and  2M)  at  a  fixed  temperature  of  85  °C.  Yang  and 
Zhao  [15]  developed  a  2D,  isothermal,  two-phase  DMFC  model. 
They  adopted  classical  multi-phase  theory  based  on  Darcy’s  law 
to  describe  two-phase  transport  in  the  anode  and  cathode  porous 
regions  and  utilized  the  drift-flux  and  homogeneous  mist-flow 
models  for  the  anode  and  cathode  flow  channels,  respectively.  In 
addition,  an  agglomerate  model  was  used  to  model  oxygen  trans¬ 
port  in  the  catalyst  layer  for  the  oxygen  reduction  reaction  (ORR). 
The  model  was  validated  against  experimental  data  that  were  col¬ 
lected  under  various  methanol  feed  concentrations  (0.25,  0.5,  1.0 
and  2.0  M).  He  et  al.  [16]  developed  a  2D,  two-phase  DMFC  model 
and  compared  the  model  predictions  with  experimental  data  for 
the  1  M  methanol  feed  concentration  and  an  operating  tempera¬ 
ture  of  70 °C.  Ge  and  Liu  [17]  developed  a  two-phase,  isothermal, 
3D  DMFC  model  and  analyzed  mainly  the  effects  of  anode  carbon 
dioxide  blockage  and  cathode  flooding.  Their  model  was  experi¬ 
mentally  validated  at  a  methanol  feed  concentration  of  1  M  and  a 
cell  temperature  of  70  °C.  Saarinen  et  al.  [18]  produced  a  3D,  single¬ 
phase  DMFC  model  and  applied  it  to  real-scale  DMFC  geometries. 
The  simulation  results  were  compared  with  current-density  dis¬ 
tribution  data  that  were  measured  with  a  segmented  cell.  In  their 
model  validation,  however,  they  adjusted  the  methanol  oxidation 
rate  at  the  cathode  catalyst  layer  to  fit  with  the  experimental  data. 


The  majority  of  the  aforementioned  DMFC  models  are  primar¬ 
ily  focused  on  methanol  crossover.  To  date,  most  of  the  DMFC 
model  validations  have  been  done  only  against  experimental  data 
that  were  under  a  narrow  range  of  operating  conditions,  particu¬ 
larly,  in  terms  of  the  methanol  feed  concentration.  Furthermore, 
few  DMFC  models  have  been  directly  compared  with  experimen¬ 
tal  data  that  are  related  to  methanol  crossover,  even  if  plenty  of 
methanol-crossover  data  are  available  in  the  literature.  Therefore, 
the  purpose  of  this  study  is  to  present  a  comprehensive  compari¬ 
son  between  numerical  simulation  results  and  experimental  data, 
using  a  ID  DMFC  model  that  was  developed  earlier  [19,20].  For 
comparison,  the  experiment  was  performed  with  reference  to  a 
unit  cell  of  25  cm2;  further,  current-density  measurements  were 
made  under  various  methanol  feed  concentrations  and  operating 
temperatures.  In  addition,  ID  simulation  results  were  compared 
with  methanol-crossover  data  that  were  measured  by  Eccarius  et 
al.  [12]  to  examine  more  precisely  methanol-crossover  phenomena 
in  DMFCs. 

2.  Experimental  details 

2.1.  MEA fabrication 

Membrane  electrode  assemblies  (MEAs)  with  an  active  area  of 
25  cm2  were  prepared  using  a  conventional  method  of  fabrica¬ 
tion  based  on  a  catalyst-coated  electrode  (CCE).  Commercial  Toray 
TGPH-060  (190  |jim  thick)  was  pretreated  with  5wt.%  PTFE  solu¬ 
tion  for  the  anode  backing  layer.  On  the  other  hand,  Sigracet  25BC 
(235  |jim  thick,  40wt.%  PTFE)  was  used  for  the  cathode  backing 
layer.  Carbon-supported  PtRu  (Hispec  12100,  Johnson  Matthey®) 
and  carbon-supported  Pt  (Hispec  13100,  Johnson  Matthey®)  were 
used  for  the  anode  and  cathode  catalysts,  respectively.  To  pre¬ 
pare  the  catalyst  slurry,  the  catalysts  were  mixed  with  5wt.% 
Nafion®  solution  (EW  1100,  Aldrich),  de-ionized  (DI)  water,  and 
isopropyl  alcohol.  After  sonication  for  10  min,  the  catalyst  slurry 
was  brushed  uniformly  on  the  backing  layers.  The  catalyst  load¬ 
ing  for  both  the  anode  and  cathode  was  2  mg  cm-2.  Nafion®  115 
membranes  (EW  1100)  from  DuPont  were  employed  on  all  the 
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Fig.  1.  Schematic  diagram  of  single-cell  DMFC  hardware. 
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Reaction  area  25cm1 2 


Fig.  2.  Schematic  diagram  of  a  three-path  serpentine  flow  channel  design  for  a 
25  cm2  test  cell  (channel  width  =  0.8  mm,  channel  depth  =  1  mm,  and  land  width  =  0.9 
mm). 


Catalyst  Coated  Membrane 


Cathode  backing  layer 


Cathode 

channel 


methanol 


water 


Anode 


Anode  backing  layer 


channel 


Fig.  3.  Schematic  of  the  computational  domain  and  species  transport  phenomena 
of  the  ID  DMFC  model. 


(3)  Perform  a  preliminary  test  three  times  in  the  constant-current 
mode  at  30,  40,  50,  60,  70  and  80  °C. 

After  the  activation  procedure,  polarization  measurements  were 
started  at  the  OCV  and  the  cell  was  operated  in  the  galvanostatic 
mode  at  a  scan  rate  of  18  mAs-1  for  each  step. 

3.  Model 


experiments.  The  membranes  were  cleaned  and  treated  by  means 
of  the  following  procedure.  The  membrane  sample  was  dipped  in 
hydrogen  peroxide  solution  (H202)  at  80  °C  for  80  min  to  eliminate 
organic  impurities  and  then  washed  with  DI  water.  Afterwards, 
it  was  immersed  in  sulfuric  acid  (H2S04)  at  80 °C  for  80 min 
to  remove  any  metallic  impurities  and  then  again  washed  with 
DI  water.  The  treated  membrane  was  positioned  between  the 
catalyst-coated  anode  and  cathode  backing  layers;  this  was  fol¬ 
lowed  by  a  hot  pressing  procedure  at  1 35  °C  at  a  pressure  of  9.8  MPa 
for  90s. 

2.2.  Cell  assembly  and  testing 

The  design  of  the  test  cell  is  shown  in  Fig.  1.  A  MEA  was 
sandwiched  between  anode  and  cathode  graphite  blocks.  Two 
gold-coated  brass  plates  were  used  as  current-collectors  for  the 
anode  and  the  cathode.  The  current-collectors  were  pressed  to 
the  graphite  blocks  by  using  end-plates  made  out  of  aluminum. 
Thermo-couples  and  electric  heaters  were  inserted  into  the  alu¬ 
minum  plates  to  maintain  the  cell  operating  temperature  at  the 
prescribed  value.  The  flow  channels  were  machined  into  the  sur¬ 
face  of  each  graphite  block  as  shown  in  Fig.  2.  The  three-pass 
serpentine  flow  channels  have  width  and  depth  of  0.8  and  1  mm, 
respectively. 

The  assembled  single  cell  was  connected  to  a  fuel  cell  test  station 
that  could  operate  in  either  a  constant-current  or  constant-voltage 
mode  and  that  was  able  to  record  current-voltage  data.  In  all 
experiments,  a  test  station  from  Fuel  Cell  Technology  was  used 
to  control  the  load,  cell  temperature,  gas  flow  rate,  humidity,  and 
back-pressure.  Methanol  solution  for  the  anode  was  supplied  using 
a  peristaltic  pump  (REGLO-CPF,  ISMATEC)  and  the  air  flow  rate  for 
the  cathode  was  adjusted  by  a  mass  flow  controller  in  the  fuel  cell 
test  station.  The  single  cell  was  activated  before  each  experiment 
as  follows. 

(1)  Flush  1  M  methanol  over  12  h  at  the  anode  side  (room  temper¬ 
ature,  2  seem),  while  the  cathode  is  closed. 

(2)  Feed  air  at  the  cathode  side  and  check  the  open-circuit  voltage 
(OCV)  for  approximately  2  h  at  70  °C. 


3.1.  Model  description  and  assumptions 

In  this  study,  a  two-phase,  isothermal,  steady-state  numerical 
DMFC  model  previously  developed  by  Ko  et  al.  [19]  is  used.  Fig.  3 
shows  the  ID  computational  domain  of  a  DMFC  MEA  that  has 
backing  layers  on  both  the  anode  and  a  cathode  sides  and  catalyst- 
coated  membrane  (CCM);  the  detailed  species  transfer  process  is 
also  given.  The  methanol  that  reaches  the  anode  catalyst  layer  (CL) 
in  the  liquid  and  the  gas  phases  dissociates  into  electrons  and  pro¬ 
tons  via  the  following  oxidation  reaction: 

MOR:  CH3OH  +  H20  ->  6H+  +  6e-  +  C02  (la) 

The  gas-phase  oxygen  passes  through  the  cathode  backing  layer 
(BL)  before  reaching  the  cathode  CL,  where  it  reduces  to  water  via 
the  following  reaction. 

ORR  :  6H+  +  6e“  +  |o2  -*  3H20  (lb) 

Due  to  the  sluggish  MOR  reaction  at  the  anode  CL  and  the 
inherent  characteristics  of  polymer  electrolyte  membranes,  such  as 
Nafion®  type  membranes  that  readily  absorb  water  and  methanol, 
methanol  crossover  is  unavoidable  in  DMFCs.  The  crossed-over 
methanol  results  in  an  additional  oxidation  reaction  at  the  cathode 
CL  that  causes  a  mixed-cathode  overpotential  and  thereby  reduces 
the  fuel  efficiency. 

The  present  model  invokes  the  following  assumptions. 

(1)  The  gas  phase  obeys  the  ideal  gas  law  because  all  the  gases  in 
a  DMFC  are  maintained  under  low  pressure  (compared  with  to 
their  respective  critical  pressures)  and  high  temperature  (above 
room  temperature). 

(2)  A  temperature  gradient  along  the  cell  thickness  is  neglected  in 
the  ID  model  due  to  the  thin  MEA  configuration  in  a  DMFC. 

(3)  An  isotropic  and  homogeneous  porous  backing  layer  (BL)  is 
assumed  and  characterized  by  effective  porosity  and  perme¬ 
ability. 

(4)  The  effect  of  C02  blockage  on  cell  performance  is  neglected  due 
to  lack  of  experimental  data  for  the  movement  of  C02  bubbles 
in  anode  BLs.  In  the  literature,  the  dynamic  behaviour  of  C02 
bubbles  has  been  mainly  observed  in  anode  channels  which 
are  not  included  in  the  computational  domain  of  the  present 
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ID  DMFC  model.  Under  this  assumption,  only  tiny  CO2  bubbles 
form  and  travel  easily  through  the  anode  BL,  which  results  in 
a  small  fraction  of  carbon  dioxide  at  the  anode  BL.  The  predic¬ 
tion  of  the  real  impact  of  C02  blockage  on  DMFC  performance 
requires  more  rigorous  studies  that  are  related  to  the  bubble 
detachment  size/frequency  and  bubble  surface  coverage,  which 
are  beyond  the  scope  of  this  work. 

(5)  Complete  consumption  of  methanol  at  the  cathode  CL  is 
assumed  after  it  crosses  over  the  membrane  from  the  anode 
to  the  cathode. 

(6)  The  electrolyte  membrane  is  fully  hydrated  due  to  the  dominant 
liquid  phase  in  DMFCs.  This  is  based  on  the  fact  that  an  aqueous 
methanol  solution  is  supplied  as  the  fuel  at  the  anode  side  and 
significant  flooding  usually  occurs  at  the  cathode  side  due  to 
the  combined  effects  of  methanol  crossover  and  the  ORR. 

(7)  Given  that  gas  permeation  through  the  electrolyte  mem¬ 
brane  is  usually  small  in  a  DMFC,  the  effects  of  oxygen 
and  carbon  dioxide  crossover  through  the  membrane  are 
neglected. 


3.2.  Species  balance  equations 


Two-phase,  species  fluxes  in  the  through-plane  direction  in  the 
anode  BL,  catalyst  membrane  (CCM)  and  cathode  BL  can  be  obtained 
by  considering  the  balances  of  individual  species,  i.e.,  methanol, 
water,  oxygen  and  carbon  dioxide  species,  which  are  expressed  in 
Eqs.  (2)  through  to  (6).  Flere,  the  species  flux  direction  from  the 
anode  to  cathode  is  defined  as  the  positive  direction. 


Methanol  balance  between  anode  BL  and  CCM: 

=  Qp  +  N M,mem 

Water  balance  between  anode  BL  and  CCM: 

Mv,a  =  Qp  +Nw,mem 

Water  balance  between  cathode  BL  and  CCM: 

^Wx  =  2f  2NM,mem  +  N\V,mem 

Carbon  dioxide  balance  between  anode  BL  and  CCM: 

Ng  —  —  — 

co2  ~  6 F 

Oxygen  balance  between  cathode  BL  and  CCM: 

No2=~4f~  2NM’mem 


(2) 

(3) 

(4) 

(5) 

(6) 


where,  I,  NM.mem .  and  NWmem  denote  the  operating  current  density, 
and  the  methanol  and  water  crossover  fluxes  through  the  mem¬ 
brane,  respectively. 


3.2.1.  Methanol  transport  equation 

The  two-phase  methanol  transport  in  the  anode  BL  is  derived 
using  Maxwell-Stefan’s  multi-component  diffusion  equation 
[21]: 


(7) 


anode  BL.  The  resultant  effective  diffusivity  can  be  written  as 
follows. 


DMff=DM^r  (8) 

Now,  Eq.  (7)  can  be  rewritten  for  the  binary  liquid-transport  system 
in  the  anode  BL  that  includes  methanol  and  water  as  follows. 


yl  jql  _  v l  jql 
KMnw,a  AWnM,a 

clDleff 

LtuM 


(9) 


In  the  above  equation,  the  molar  fractions  of  methanol  and  water 
are  defined  as  follows. 


— 


_M  . 
rl  ’ 


rl 


-  c_w_  _  1  _  xi 

r,  w  -  l  -  1  XM 

4  4 


(10) 


where  c\  represents  the  total  molar  concentration  of  the  liquid 
phase  in  the  anode  BL,  i.e.,  the  sum  of  the  individual  molar  con¬ 
centrations,  namely: 

Ct  —  CM  +  CW  (11) 


Using  Eqs.  (10)  and  (11),  Eq.  (9)  for  the  methanol  flux  through 
the  anode  BL,  Nj^  fl,  can  be  rewritten  as  follows. 


Nl  — 


rl  Ml  -  ClDl,e^ VC1 
LMiVW,a  LtuM  VLM 


4  -  C 


M 


(12) 


The  gas-phase  transport  in  the  anode  BL  is  assumed  to  be  a 
binary  diffusion  system  that  consists  of  methanol  vapour  and  car¬ 
bon  dioxide.  Again,  using  the  Maxwell-Stefan  equation,  viz.,  Eq.  (9), 
the  molar  fraction  of  methanol  vapour  can  be  written  as: 


Y7vS  _  M  CO 2, a  ACO 2  M, a 

4 


(13) 


In  Eq.  (13),  the  effective  diffusivity  of  methanol  vapour, 
in  the  anode  BL  can  be  a  function  of  the  porosity,  liquid  saturation, 
and  the  Bruggeman  factor,  n,  as  follows. 


D^  =  f4[£(l-s)]" 


(14) 


According  to  the  ideal  gas  law,  the  total  molar  concentration  of 
gas,  cf ,  can  be  determined  as  follows. 


(15) 


Now,  Eq.  (13)  can  be  written  for  the  methanol  vapour  flux  along 
the  anode  BL,  Nj^  fl,  as: 


Ng  -  r8 D8,6^ Vc8 
mS  _  M  CO  2, a  4  UM  VCM 

iVM,a  -  rg 


(16) 


Assuming  methanol  and  methanol  vapour  to  be  in  thermodynamic 
equilibrium  and  applying  Henry’s  law  to  Eq.  (16),  we  can  express 
the  methanol  vapour  flux,  N8^  fl,  in  terms  of  clM  instead  of  c 8/[: 


N8  - 

iVM,a  — 


rl  N8  -c^D^’^Vc* 
4rvC02,a  LtuM  VCM 


kHcf  -  cj 


M 


(17) 


where  Henry’s  constant,  kH,  is  defined  as  the  ratio  between 
the  methanol  concentration  in  the  liquid  and  gas  phases  and  is 
expressed  as  a  function  of  temperature,  i.e., 


For  a  binary  liquid-transport  system  through  the  porous  BLs 
of  DMFCs,  the  binary  diffusivity,  Dap,  should  be  modified  via  the 
Bruggeman  correlation  to  account  for  the  effect  of  the  porosity, 
tortuosity  of  the  anode  BL,  and  liquid  saturation,  i.e.,  the  vol¬ 
ume  fraction  of  liquid  that  occupies  the  open  pore  space  of  the 


kH  =  CJ  =  k0e-(AHso,/R)((l/T)-l/T0)Rr  (18) 

4w 

where,  k ^  refers  to  Henry’s  constant  at  the  reference  state  of 
T0  =  298.15  K  and  AHso(  refers  to  the  enthalpy  of  the  solution. 
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Using  Eqs.  (2),  (3),  (5),  (12),  and  (17),  it  is  possible  to  obtain  the 
final  form  of  the  methanol  balance  equation  between  the  anode  BL 
and  the  membrane  as  follows: 


c^((f/6F)  +  Nw,mem)-c(DjMv4  |  +(-J/6F)-cfl4v+ 

Ct  ~  CM  kf/Cf  -  + 

—  7TF  +  hi m, mem  (19) 


In  Eq.  (19),  the  methanol  flux  that  passes  through  the  membrane, 
N M,mem,  can  be  driven  by  the  EOD  due  to  the  proton  flux  and  dif¬ 
fusion  that  arise  from  the  methanol  concentration  gradient  across 
the  membrane: 


Nm, 


,mem  ^M,mem  +  ^M,mem  ^d,M  p  +  Dm,i 


I  aCL 


(20) 


where  nd  M  denotes  the  EOD  coefficient  of  methanol  given  by  Liu 
[22],  which  is  a  function  of  the  EOD  coefficient  of  water,  ndw,  and 
the  methanol  and  water  concentrations  in  the  anode  catalyst  layer. 
That  is, 


nd,M  =  nd,w  -r- 
cw 


(21) 


32.2.  Oxygen-transport  equation 

The  diffusion  flux  of  oxygen  through  the  cathode  BL  is  expressed 
using  Fields  diffusion  equation.  Therefore,  with  Eqs.  (6)  and  (20),  the 
final  form  of  the  oxygen  transport  equation  for  the  cathode  side  can 
be  obtained  as  follows: 

NS02  =  =  -i  - 

=  -(6 nd)M  +  1)T  _  ^DM  mem^LcL  (22) 

4L  2  Ojnem 

where  is  the  effective  diffusion  coefficient  of  oxygen  through 
the  cathode  BL  and  is  expressed  as  a  function  of  the  porosity,  liquid 
saturation,  and  Bruggeman  factor,  n,  as  follows. 


hydrophilic  porous  layers  as  follows. 

J  1.417(1  -  s)- 2.120(1  -s)2  +  1.263(1  -s)3  if#  <  90° 

1  1.417s- 2.120s2  +  1.263s3  if#  >  90° 

(27) 


The  liquid-phase  Darcy’s  equation,  Eq.  (25),  can  be  rewritten  in 
terms  of  the  capillary  pressure,  Pc,  by  assuming  that  the  pressure 
drop  through  the  gas  flow  is  relatively  negligible  when  compared 
with  the  gradient  of  the  liquid  pressure. 


,J_  K,Jv7d  _  K.,,  2a  cos  OdJ 

ds 


Ul=  -r/4VPc  =-X- 


Vs 


(28) 


By  combining  Eqs.  (12)  and  (28)  into  Eq.  (24),  the  final  form  of  the 
liquid-transport  equation  at  the  anode  can  be  expressed  as  follows: 


,  I<a  ,  /  2a  •  cos  #  dj  ^ 

P—A - - - ^Vs  =  -Mm 


rl  Nl  -  rl  D‘ 
LMnW,a  LtuM 


ds 


cl  -cl 

ct 


—  MwNlWa  (29) 

where  Nlw  a  denotes  the  water  flux  through  the  anode  BL.  According 
to  Eq.  (3),  Nlw  a  is  equivalent  to  the  sum  of  the  water  consumption 
rates  through  the  MOR  at  the  anode  side,  //6F,  and  the  water  flux 
across  the  membrane,  NWjiriem,  i.e.,  driven  by  the  EOD  due  to  the 
proton  flux,  the  diffusion  due  to  the  water-content  gradient,  and 
the  hydraulic  pressure  gradient,  as  shown  by: 

]\r.A/  -  mEOD  ,  Ndiff  _  Npl 

iV W,mem  -  ^W,mem  +  lSW,mem  nW,mem 

_  I  n  Pmem  ha  —  hc  Kmem  P£  —  P a 

—  Eld,W  +  +  El w,  mem  “+777  e  71  J 

F  EW  dmem  MWVl  Omem 

Similarly,  we  can  derive  the  liquid-transport  equation  at  the  cath¬ 
ode  side  as  follows. 


jKClJ  2a  cos  0dj(s) 

pJJkr~  dT 


2p  T  2NM,mem 


+  %,n 


(31) 


Dgf  =  DS2[e(l-s)]n  (23) 

32.3.  Liquid-saturation  equation 

The  derivation  of  the  liquid  saturation  equation  commences 
from  the  mass  conservation  equation  for  the  liquid  phase.  For  the 
anode  side,  the  mass  conservation  equation  for  the  liquid  phase  can 
be  expressed  as: 

p'u1  =  -MMN'M  a  -  MwNlW  a  (24) 

On  the  other  hand,  Darcy’s  law  can  be  employed  to  describe 
liquid  flow  in  porous  media,  such  as  the  anode  BL: 


3.2.4.  Electrochemical  kinetics 

At  a  given  current  density,  /,  the  cell  voltage  is  determined  from 
the  theoretical  thermodynamic  potential,  anode  and  cathode  over¬ 
potentials,  membrane  ohmic  resistance,  and  contact  resistance  as 
follows. 

Vcell  =  U° -ria  +  ric-IRcont-I—  (32) 

Kmem 

Regarding  Eq.  (32),  it  should  be  mentioned  that  the  theoretical  equi¬ 
librium  potential  U°  represents  the  difference  between  the  anode 
and  cathode  standard  potentials  (U°  =  U°  -  U£),  which  are  respec¬ 
tively  determined  by  using  the  Nernst  equation  as  follows. 


ul  =  -F-k'VP1  (25) 

where  I<  and  k[  denote  the  permeability  of  the  given  porous 
medium  and  the  relative  permeability  for  the  liquid  phase, 
respectively.  The  capillary  pressure,  Pc  can  be  expressed  as 
a  function  of  the  surface  tension,  o ,  porosity,  £,  permeabil¬ 
ity,  I<,  and  contact  angle,  #,  of  the  porous  medium  as  follows 
[23]: 

Pc  =  ps -P1  =acos  0(~)  J(s)  (26) 

where  the  Leverett  function,  J(s),  is  the  dimensionless  capil¬ 
lary  pressure,  which  can  be  expressed  for  both  hydrophobic  and 


For  the  anode: 

—  A  2"® 

U°a  =  „-fa  =  -0.016  V  vs.  standard  hydrogen  electrode  (SHE) 

br 

(33) 

For  the  cathode: 

-As0 

U°  =  =  1.229V  vs.  SHE  (34) 

6  r 

The  second  and  third  terms  in  the  right-hand  side  of  Eq.  (32)  are 
the  anode  kinetic  loss,  qa,  and  the  cathode  kinetic  loss,  rjc-  These 
two  parameters  are  calculated  from  the  electrochemical  kinetics 
expressions  obtained  from  the  Butler-Volmer  equation. 
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For  the  anode  MOR: 

,  =  -ref  CM\acLexp((YalRT)Fr}a)  . 

0’aCM\aCL+KreacexpttYalRT)Fr)a )  1  J 

For  the  cathode  ORR: 

/  + Wr  =  io^d  -  Sc)  j  exp  (H^c)  (36) 

The  MOR  kinetic  expression  in  Eq.  (35)  was  derived  by  Mey¬ 
ers  and  Newman  [24]  from  a  four-step  mechanism  proposed  by 
Gasteiger  et  al.  [25].  In  Eq.  (35),  cM\acL  and  I<reac  stand  for  the 
methanol  concentration  at  the  anode  CL  and  the  reaction  constant, 
respectively.  On  the  other  hand,  the  ORR  kinetics  expression  in  Eq. 
(36)  is  based  on  first-order  Tafel  kinetics  and  is  obtained  from  the 
experimental  studies  of  Bernardi  and  Verbrugge  [26],  and  Gottes- 
feld  and  Zawodzinski  [27].  In  addition,  ir^a  and  denote  the 
exchange  current  densities  for  the  anode  MOR  and  cathode  ORR, 
respectively,  which  are  given  by  Wang  and  Wang  [13]  as  follows. 


the  methanol  and  oxygen  concentrations  in  the  anode  and  cathode 
channel  inlets,  respectively. 


Methanol  concentration  for  anode: 


Oxygen  concentration  for  cathode: 


In  the  above,  §a  and  §c  represent  the  anode  and  cathode  stoi¬ 
chiometric  numbers,  respectively. 

On  the  other  hand,  for  the  liquid-transport  equations,  viz.,  Eqs. 
(29)  and  (31),  the  values  of  the  interfacial  liquid  coverage  of  the 
anode  and  the  cathode  BL  surfaces  are  assumed  to  be  0.8  and  0.1, 
respectively. 

For  the  anode: 


For  the  anode  MOR: 


For  the  cathode  ORR: 
ir0efc{T)  =  0.0422  exp 


f 35570  ( 

(  1 

--)) 

1  R  ' 

v353 

tJJ 

( 73200 

fJL 

'  V  R 

V353 

Tj) 

(37) 


(38) 


Finally,  the  last  two  terms  in  Eq.  (32)  are  related  to  the  ohmic 
losses  due  to  the  contact  resistance  and  membrane  resistance, 
respectively.  Regarding  the  membrane  resistance,  8mem  and  Kmem 
stand  for  the  membrane  thickness  and  proton  conductivity  through 
the  membrane,  respectively. 


3.3.  Boundary  conditions 

The  species  concentrations  that  are  applied  to  the 
anode | cathode  BL| channel  interfaces  are  calculated  based  on 


S\aBC  =  0-8 

(41) 

For  the  cathode: 

SlcBC  =  0-1 

(42) 

The  above  assumptions  are  based  on  the  fact  that  in  DMFCs, 
the  two-phase  flow  in  the  anode  channel  is  more  like  a  pure  liquid 
flow  with  a  small  gas  fraction,  whereas  the  gas  phase  is  usually 
more  dominant  in  the  two-phase  flow  pattern  along  the  cathode 
gas  channel. 

3.4.  Numerical  procedures 

The  governing  equations  derived  in  the  foregoing  section  are 
implemented  using  a  FORTRAN  code  and  solved  simultaneously 
in  the  entire  computational  domain  that  comprises  the  anode  BL, 
CCM,  and  cathode  BL.  After  the  computational  domain  is  initialized, 
the  methanol  transport  equation  is  first  solved.  Then,  the  cell  volt¬ 
age  is  calculated  from  the  known  current  density  and  methanol 


Table  1 

Transport  and  physical  properties. 


Description 

Value 

Ref. 

Water  density,  pw 

1.0  x  103  kgm-3 

[22] 

Methanol  density,  pM 

0.79  x  103  kgm-3 

[22] 

Membrane  density,  pmem 

1.98  x  103  kgm-3 

[22] 

Faraday  constant,  F 

96,  487  C  mol-1 

[13] 

Membrane  equivalent  weight,  EW 

1.1  kg  mol-1 

[22] 

Porosity  of  anode  backing  layer,  ea 

0.6 

[30] 

Porosity  of  cathode  backing  layer,  ec 

0.6 

[30] 

Cell  contact  resistance,  Rcont 

0.2  x  10-4  £2  m2 

Assumed 

Anode  BL  permeability,  Kqbl 

1.0  x  10-12  m2 

[15] 

Cathode  BL  permeability,  I<cbl 

1.0  x  10-12  m2 

[15] 

Membrane  permeability,  /Cmem 

4.0  x  10-20  m2 

[30] 

Surface  tension,  a 

0.0625  Nm-1 

[16] 

Proton  conductivity  of  the  membrane,  Kmem 

(0.5139A.- 0.326)  exp  [l268  (^  - 

0] 

sm_1 

[13] 

"  0.458509  -  5.34047  x  10“3T+ 

Water  viscosity,  p} 

2.31231  x  10“5T2  -4.49161  x  10 

— 8  y\3 

kg  m_1  s_1 

[13] 

_ +3.27681  x  10“nT4 

Methanol  diffusivity  in  liquid,  DlM 

IQ-5.4163-999.778/T  m2  S~1 

[13] 

Methanol  diffusivity  in  gas, 

-6.954  x  10“2  +4.5986  x  10“4T 

x  10 

r4  m2  s_1 

[13] 

+9.4979  x  10“7T2 

Oxygen  diffusivity  in  gas, 

—  /  \  1.823 

1.775  x  10~5(  273  J  m2  s-1 

[13] 

Methanol  diffusivity  in  the  membrane,  DM,mem 

1.98  xlO-10  exp  (2416  (2k  -  l)] 

|  m2s 

-i 

Fitted  [31] 

Water  diffusivity  in  the  membrane,  Dw,mem 

4.80  x  10“n  exp  (2416  (3J3  -  7)) 

m2  s 

-i 

[30] 

Electro-osmotic  drag  coefficient  of  water,  ndiW 

1.6767  +  0.0155  x  (T- 273)  +  8.9074  ; 

<(T- 

273)2 

Fitted  [32] 
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crossover  rate.  Next,  the  liquid-saturation  and  oxygen-transport 
equations  are  solved  successively.  After  all  the  transport  equations 
are  calculated,  the  related  parameters  are  updated.  The  iterations 
proceed  until  the  methanol  concentration  at  the  anode  CL  satisfies 
the  convergence  criterion  (<10-5). 

4.  Results  and  discussion 

For  the  present  model  validation  study,  the  polarization  curve 
measurements  were  conducted  using  an  in-house  25  cm2  MEAs. 
The  voltage-current  density  data  were  collected  at  three  different 
cell  temperatures  (50,  60,  and  70  °C)  and  four  different  methanol 
feed  concentrations  (1,  2,  3  and  4M).  For  all  experimental  cases, 
the  flow  rates  for  the  anode  and  cathode  sides  were  fixed  at 
8  seem  for  the  methanol  solution  and  1500  seem  for  air;  these 
lead  to  sufficiently  high  stoichiometric  factors  for  the  anode  and 
cathode  flow  rates  (roughly  30  for  both  the  anode  and  cathode 
sides  based  on  1M  methanol  and  dry  air  at  the  reference  cur¬ 
rent  density  of  100mA cm-2).  It  should  be  noted  that  the  use 
of  a  sufficiently  high  air  flow  rate  in  the  present  experiment  is 
to  prevent  the  generation  of  negative  (electrolytic)  current,  i.e., 
caused  by  methanol  crossover  followed  by  MOR  at  the  cathode  side 
[28,29]. 

For  a  comparison  of  the  results  between  the  simulations  and 
experiments,  one-dimensional  numerical  simulations  were  carried 
out  under  the  same  geometric  (through-plane  direction)  and  oper¬ 
ating  conditions  as  in  the  experiments.  It  should  be  emphasized 
that  all  the  input  parameters  listed  in  Table  1  were  used  for  all  the 
simulation  cases  and  no  adjustable  parameters  were  considered  for 
fitting  with  the  experimental  data. 

A  comparison  of  the  calculated  and  measured  polarization 
curves  at  three  operating  temperatures  (50,  60  and  70  °C)  for  a 
methanol  feed  concentration  of  1  M  is  given  in  Fig.  4.  Gener¬ 
ally,  good  agreement  between  the  simulations  and  experiments  is 
achieved  except  in  the  region  of  low  current  density.  The  deviations 
can  be  attributed  to  the  use  of  the  Tafel  kinetic  expressions,  viz., 
Eqs.  (35)  and  (36),  and  the  kinetic  parameters,  namely,  Eqs.  (37)  and 
(38).  The  effect  of  the  operating  temperature  is  seen  clearly  in  Fig.  4. 
Both  the  calculated  and  the  measured  cell  performance  increase 
with  the  temperature  within  the  50-70  °C  range.  This  implies  that 
the  positive  effect  due  to  the  improved  electrochemical  kinetics  and 
reactant  transport  with  rising  temperature  is  more  prominent  than 
the  negative  effect  caused  by  methanol  crossover  that  is  enhanced 
by  the  rising  temperature. 


Fig.  4.  Comparison  of  the  simulated  (lines)  and  measured  (symbols)  polarization 
curves  at  different  cell  temperatures  (50, 60,  and  70  °C)  for  a  methanol  concentration 
of  1  M. 


Fig.  5.  Comparison  of  the  simulated  (lines)  and  measured  (symbols)  polarization 
curves  at  different  cell  temperatures  (50, 60,  and  70  °C)  for  a  methanol  concentration 
of  2  M. 

Simulated  and  measured  polarization  curves  for  a  methanol 
feed  concentration  of  2  M  are  presented  in  Fig.  5.  Again,  excellent 
agreement  between  the  results  of  the  simulations  and  those  of  the 
experiments  is  seen,  although  there  is  inadequate  agreement  at 
low  current  densities  due  to  the  Tafel  approximation.  When  com¬ 
pared  with  the  case  of  1  M  (Fig.  4),  the  polarization  curves  for  2  M 
(Fig.  5)  reveal  a  reduction  in  the  mass-transport  limitation  due  to 
the  use  of  a  higher  methanol  feed  concentration.  Therefore,  the 
limiting  current  densities  in  the  cases  of  2  M  exist  in  the  vicinity  of 
0.45-0.5  A  cm-2,  whereas  those  for  1  M  cases  are  0.2-0.3  A  cm-2. 

Simulated  polarization  curves  with  the  experimental  data  for  a 
methanol  feed  concentration  of  3  M  are  given  in  Fig.  6.  Similar  to 
Figs.  4  and  5,  good  agreement  between  the  results  of  the  simula¬ 
tions  and  those  of  the  experiments  is  observed.  When  compared 
with  Fig.  5,  the  data  in  Fig.  6  clearly  show  that  the  cell  perfor¬ 
mance  for  3  M  is  slightly  lower  than  for  2  M,  which  is  indicative 
of  the  enhanced  methanol  crossover  due  to  the  raised  methanol 
feed  concentration.  The  results  of  model  validation  for  a  methanol 
feed  concentration  of  4  M  are  given  in  Fig.  7.  Once  again,  similar 
to  low  methanol  feed  concentrations,  by  comparing  the  results 
of  the  simulations  with  those  of  the  experiments,  it  is  seen  that 
the  model  validation  is  satisfactory.  Nevertheless,  the  agreement  is 
relatively  poor  at  70  °C  where  the  model  overestimates  cell  per¬ 
formance.  It  should  be  noted  that  the  differences  between  the 


Fig.  6.  Comparison  of  the  simulated  (lines)  and  measured  (symbols)  polarization 
curves  at  different  cell  temperatures  (50, 60,  and  70  °C)  for  a  methanol  concentration 
of  3  M. 


J.  Ko  et  al.  /  Journal  of  Power  Sources  196(2011)  935-945 


943 


Fig.  7.  Comparison  of  the  simulated  (lines)  and  measured  (symbols)  polarization 
curves  at  different  cell  temperatures  (50, 60,  and  70  °C)  fora  methanol  concentration 
of  4M. 


polarization  curves  at  the  three  different  cell  temperatures  are 
reduced  more  in  the  case  of  4  M  than  in  the  other  cases  (of  1  -3  M). 
This  implies  that  for  4M,  the  negative  effect  due  to  enhanced 
methanol  crossover  by  the  rising  temperature  becomes  comparable 
with  the  positive  effect  that  is  due  to  the  enhanced  electrochemi¬ 
cal  reaction  kinetics  and  reactant  transport  with  rising  temperature 
(Fig.  7). 

The  present  ID  DMFC  model  was  also  validated  against  the 
methanol  crossover  current  density  data  that  were  measured  by 
Eccarius  et  al.  [12].  Simulated  and  measured  crossover  current 
density  curves  as  a  function  of  the  operating  current  density  is  pre¬ 
sented  Fig.  8.  It  should  be  noted  that  the  crossover  current  density 
can  be  calculated  according  to  Faraday’s  law  via  Eq.  (20)  as  follows. 

cl  I 

fxover  —  6FNm 

,mem  —  6nd^I  +  6FDm  ,mem  (43) 

v _ _ /  Omem 

s _ _ _ ✓ 

EOD  v 

diffusion 

The  crossover  current  density  is  decomposed  into  the  EOD 
and  diffusion  components.  The  conditions  for  the  simulations  and 
experiments  are  based  on  a  methanol  feed  concentration  of  1 .5  M,  a 


Fig.  8.  Simulated  (lines)  and  measured  (symbols)  methanol  crossover  current  den¬ 
sity  curves  as  a  function  of  the  cell  current  density  and  the  contributions  of  the 
diffusion  and  electro-osmotic  drag  (EOD)  to  the  methanol  crossover  current  density 
(Nation  117  with  a  catalyst  loading  of  2.5  mg  cm-2  on  both  electrodes,  50  °C,  1.5  M 
solution  of  methanol,  anode  flow  rate  of  3  ml  min-1 ,  and  the  cathode  stoichiometry 
of  6). 


Fig.  9.  Simulated  (lines)  and  measured  (symbols)  methanol  crossover  current  den¬ 
sity  curves  as  a  function  of  the  cell  current  density  and  the  contributions  of  the 
diffusion  and  electro-osmotic  drag  (EOD)  to  the  methanol  crossover  current  density 
(Nation  1135  with  a  catalyst  loading  of  2.5  mg  cm-2  on  both  electrodes,  0.5  M  solu¬ 
tion  of  methanol,  anode  flow  rate  of  1  ml  min-1 ,  and  the  cathode  stoichiometry  of  2 
for  different  operating  temperatures). 


cell  temperature  of  50  °C,  and  a  Nation®  117  membrane  ( 1 80  pun). 
There  is  good  agreement  between  the  experimental  data  (symbols) 
and  simulation  results  (lines);  both  show  an  increasing  trend  with 
the  cell  current  density.  The  simulation  results  further  present  the 
contribution  of  the  diffusion  and  EOD  to  the  total  crossover  cur¬ 
rent  density.  At  the  OCV,  the  effect  of  the  EOD  is  not  present  and 
thus,  the  methanol  crossover  is  fully  driven  by  diffusion.  As  the 
operating  current  density  increases,  the  effect  of  diffusion  con¬ 
tinuously  decreases  because  the  methanol  concentration  in  the 
anode  catalyst  layer  is  lowered  due  to  the  increasing  rate  of  the 
MOR.  On  the  other  hand,  according  to  Eq.  (43),  the  EOD  component 
of  the  methanol  crossover  current  density  continues  to  increase 
with  cell  current  density.  Therefore,  the  overall  increase  in  the 
total  crossover  current  density  with  increasing  cell  current  density 
implies  that  the  increasing  rate  of  the  EOD  component  is  more  dom¬ 
inant  than  the  decreasing  rate  of  the  diffusion  component  under  the 
given  conditions. 

A  comparison  of  the  simulated  and  measured  methanol 
crossover  current  densities  at  the  two  operating  temperatures  of 
30  and  50  °C  is  seen  in  Fig.  9.  One  more  case  corresponding  to  a  cell 
temperature  of  70  °C  was  additionally  simulated  to  investigate  the 
influence  of  the  operating  temperature  on  the  methanol-crossover 
phenomena.  All  these  cases  are  based  on  a  methanol  feed  con¬ 
centration  of  0.5  M  and  a  Nation®  1135  membrane  (90p,m).  The 
curves  of  methanol  crossover  current  density  match  well  across 
the  simulations  and  experiments  and  that  both  the  measured 
and  simulated  crossover  current  densities  decrease  the  cell  cur¬ 
rent  density,  which  indicates  that  the  crossover  current  densities 
under  these  conditions  are  mainly  determined  by  diffusion.  In 
addition,  as  the  cell  temperature  is  elevated,  the  diffusion  compo¬ 
nent  of  the  crossover  current  density  increases  significantly,  which 
illustrates  that  the  diffusion  component  is  the  main  contributor 
to  the  rise  in  the  crossover  current  density  with  the  tempera¬ 
ture.  According  to  the  methanol  crossover  diffusivity  correlation, 
Dm, mem  in  Table  1,  it  is  obvious  that  the  methanol  at  the  anode 
can  more  easily  diffuse  through  the  membrane  at  higher  cell 
temperatures. 

It  is  expected  that  the  rate  of  methanol  crossover  can  be  affected 
by  the  thickness  of  the  membrane.  The  crossover  current  den¬ 
sity  curves  simulated  with  three  different  membranes,  i.e.,  Nation® 
117  (180  p,m),  Nation®  1135  (90  pan),  and  Nation®  112  (50  pan), 
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Fig.  10.  Simulated  (lines)  and  measured  (symbols)  methanol  crossover  current  den¬ 
sity  curves  as  a  function  of  the  cell  current  density  and  the  contributions  of  the 
diffusion  and  electro-osmotic  drag  (EOD)  to  the  methanol  crossover  current  den¬ 
sity  (a  catalyst  loading  of  2.5  mg  cm-2  on  both  electrodes,  50  °C,  0.5  M  solution  of 
methanol,  anode  flow  rate  of  1  ml  min-1,  and  the  cathode  stoichiometry  of  6  for 
different  membrane  thicknesses). 


are  presented  in  Fig.  10.  Nation®  117  and  Nation®  1135,  the  pre¬ 
dicted  results  can  be  compared  with  the  experimental  data.  Again, 
these  crossover  current  density  curves  demonstrate  good  agree¬ 
ment  between  the  results  of  the  simulations  and  those  of  the 
experiments.  Also,  the  results  of  the  simulation  model  further  indi¬ 
cate  that  the  diffusion  component  of  the  crossover  current  density 
is  significantly  reduced  with  thicker  membranes,  whereas  the  EOD 
component  remains  almost  constant,  regardless  of  the  membrane 
thickness. 

The  final  results  of  validation  are  presented  in  Fig.  1 1  for  0.5 
and  1 .5  M  methanol  feed  concentrations.  Once  again,  there  is  excel¬ 
lent  agreement  in  the  crossover  current  density  curves  between 
the  simulations  and  experiments  for  both  0.5  and  1.5  M,  and 
the  model  successfully  captures  the  experimental  trends  in  the 
crossover  current  density  under  various  methanol  feed  concen¬ 
trations.  In  particular,  it  should  be  noted  that  a  change  in  the 
methanol  feed  concentration  leads  to  a  significant  change  in  the 
shape  of  the  crossover  current  density  curve,  as  well  as  in  the 
degree  of  methanol  crossover.  For  the  0.5  M  case,  the  crossover 
current  densities  are  relatively  low  and  decrease  with  cell  cur¬ 
rent  density,  whereas  the  1.5  M  case  exhibits  higher  levels  of 
methanol  crossover  and  an  increasing  trend  in  the  crossover  cur¬ 
rent  density  with  the  cell  current  density.  The  detailed  simulation 
results  demonstrate  that  both  the  EOD  and  diffusion  crossover 
currents  increase  significantly  as  the  methanol  feed  concentration 
increases.  In  particular,  the  EOD  crossover  current  density  increases 
dramatically  with  cell  current  density  under  high  methanol  feed 
concentrations  which,  as  seen  in  Eq.  (43),  results  from  increases 
in  both  the  methanol  EOD  coefficient,  ndM  and  the  cell  current 
density,  I.  For  the  successful  use  of  methanol  fuel  with  high  con¬ 
centration  for  DMFCs,  it  is  well  known  that  the  development  of 
membranes  that  have  greater  resistance  to  methanol  crossover  is 
important  because  the  diffusion  component  of  methanol  crossover 
can  be  significantly  reduced  by  such  membranes.  Nevertheless, 
the  simulation  results  in  Fig.  1 1  indicate  further  that  the  reduc¬ 
tion  of  the  methanol  concentration  in  the  anode  catalyst  layer  is 
even  more  important  in  order  to  suppress  simultaneously  both 
the  diffusion  and  EOD  components.  This  reduction  can  be  achieved 
by  designing  a  highly  resistive  anode  backing  layer  for  methanol 
transport. 


Fig.  11.  Simulated  (lines)  and  measured  (symbols)  methanol  crossover  current  den¬ 
sity  curves  as  a  function  of  the  cell  current  density  and  the  contributions  of  the 
diffusion  and  electro-osmotic  drag  (EOD)  to  the  methanol  crossover  current  density 
(Nation  117  with  a  catalyst  loading  of  2.5  mg  cm-2  on  both  electrodes,  50  °C,  anode 
flow  rate  of  1  ml  min-1,  and  the  cathode  stoichiometry  of  6  for  different  methanol 
feed  concentrations). 


5.  Conclusions 

The  main  objective  of  this  study  is  to  compare  the  simulation 
results  of  a  ID  DMFC  model  with  experimental  data  that  are  mea¬ 
sured  under  a  wide  range  of  cell  designs  and  operating  conditions.  It 
should  be  noted  that  no  fitting  parameters  have  been  assumed  for 
the  comparison  activities  and  that  the  same  model  input  param¬ 
eters  have  been  employed  for  all  the  cases  of  comparison.  This 
study  mainly  involves  a  comparison  of  the  simulated  and  mea¬ 
sured  values  for  the  cell  current  density  data,  which  were  obtained 
from  in-house  fuel  cell  tests  and  the  methanol  crossover  current 
density  data  given  by  Eccarius  et  al.  [12].  In  general,  all  the  com¬ 
parisons  yield  excellent  agreement  under  different  methanol  feed 
concentrations  (1-4M),  cell  temperatures  (50-70  °C)  and  mem¬ 
brane  thicknesses  (90  and  180  pan),  which  proves  the  validity  and 
accuracy  of  the  present  ID  DMFC  model.  In  addition,  a  detailed 
examination  is  undertaken  of  the  simulation  results  and  experi¬ 
mental  data.  The  trends  observed  in  this  study  can  be  summarized 
as  follows. 

Both  the  simulation  results  and  experimental  data  demonstrate 
clearly  that  the  operation  of  a  DMFC  at  a  high  cell  temperature  helps 
to  improve  cell  performance  due  to  the  enhanced  electrochemical 
kinetics  and  the  transportation  ability  of  reactants.  The  beneficial 
effect  of  high-temperature  operation  is  diminished,  however,  when 
highly  concentrated  methanol  fuel  is  employed  for  DMFC  opera¬ 
tions.  The  trend  is  indicative  of  enhanced  methanol  crossover  flux 
through  the  membrane  under  rising  cell  temperature. 

In  the  methanol  crossover  current  density  profiles,  the  decrease 
in  methanol  crossover  density  with  cell  current  density  can  be 
ascribed  to  the  dominance  of  the  methanol  diffusion  effect  over 
the  EOD  effect.  On  the  other  hand,  the  methanol  crossover  den¬ 
sity  increases  with  cell  current  density  when  the  EOD  effect  is 
more  dominant.  Simulations  indicate  that  the  methanol  crossover 
flux  through  diffusion  is  sensitive  to  all  the  key  DMFC  design  and 
operating  parameters  considered  in  this  study.  Methanol  crossover 
by  diffusion  is  enhanced  under  greater  methanol  feed  concen¬ 
trations  and  cell  temperatures  but  suppressed  under  a  rising 
membrane  thickness.  Unlike  methanol  diffusion  across  the  mem¬ 
brane,  methanol  crossover  that  is  driven  by  the  EOD  is  not  altered 
much  with  variation  in  cell  temperature  and  membrane  thick- 
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ness.  By  contrast,  the  methanol  feed  concentration  is  found  to  be 
a  significant  factor  and  considerably  affects  the  EOD  component 
of  methanol  crossover.  Therefore,  the  presented  simulation  results 
further  indicate  that  for  the  use  of  highly  concentrated  methanol 
fuel  to  increase  the  specific  energy  of  a  DMFC  system,  both  the 
EOD  and  diffusion  components  of  methanol  crossover  should  be 
suppressed,  which  can  be  achieved  by  designing  a  highly  resistive 
anode  backing  layer  for  methanol  transport. 
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